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Abstract. Atomic detection by fluorescence may fail because of reflection from the 
laser or transmission without excitation. The detection probability for a given velocity 
range may be improved by controlling the detuning and the spatial dependence of the 
laser intensity. A simple optimization method is discussed and exemplified. 
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Figure 1. Detection probability versus velocity v for laser beam width L — 10 /im. 
(This and other figures are obtained for the transition at 852 nm of Cs atoms, with 
7 = 33.3 x 10 6 s" 1 .) fl = 0.1033 x 10 6 s" 1 (solid line); fl = 5 7 (dashed line); in both 
cases A = 0. 

1. Introduction 

One of the standard ways to measure the time of flight, or simply the presence of an 
atom, consists of illuminating it with a laser and detect the induced fluorescence. There 
are many different experimental settings depending on the incident atomic velocities 
and spatial span of the atomic cloud or beam. In particular, the light of the probe laser 
may be spread as a broad sheet perpendicular to the atomic motion or be focused onto a 
diameter of a few microns; it may also be continuous or pulsed. These detection schemes 
may fail however at low (ultracold) atomic velocities because of atomic reflection from a 
strong laser field; the atoms may also traverse the finite laser-illuminated region without 
emitting any photon even at moderate velocities if the laser intensity is weak. 

Figure 1 illustrates these two phenomena. It represents the detection probability 
versus (ultracold) atomic velocity for a given laser-beam width L and two different 
laser intensities. This and other figures below are for Cs atoms, and the transition 
5*1/2) F = 4 — > P3/2, F — 5, with Einstein coefficient 7 = 33.3 x I0~ 6 s _1 . The minimum 
velocity considered is 0.2 cm/s, slightly below the recoil velocity ~ 0.35 cm/s. In each 
case the laser intensity has been assumed to be constant in the illuminated region for 
simplicity. Note the significant atomic reflection for strong driving and low velocities. 
By contrast the weak laser detects the very slow atoms but fails for the faster ones. 
No intermediate intensity between the two cases depicted can avoid the two problems 
simultaneously in an arbitrarily broad velocity range. An ideal detection should be able 
to avoid both effects at least for the velocity spread of the atomic clouds of interest. 
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The objective of this paper is to show that an appropriate adjustment of detuning and 
of the position dependence of the laser intensity may indeed achieve this goal. 

In the next section we shall present the theory to model the atom detection, section 
III describes the optimization of laser intensity and detuning, and in section IV we 
provide a numerical example. 



2. Basic theory 

The fundamental theory for the modeled experiment is described in Refs. pQ, [2] and 
jsJj. A two- level atom with transition frequency u impinges along the x direction on a 
perpendicular continuous laser beam of frequency ujl directed along the y direction. In 
the so called quantum jump approach [3] the continuous measurement of the fluorescence 
photons is simulated by a periodic projection onto no-photon or one-photon subspaces 
every At, a time interval large enough to avoid the Zeno effect, but smaller than any 
other characteristic time. The amplitude for the wavepacket of undetected atoms in 
the interaction picture for the internal Hamiltonian obeys, in a time scale coarser than 
At, and using the rotating wave and dipole approximations, an effective Schrodinger 
equation governed by the complex "conditional" Hamiltonian (the hat is used to 
distinguish momentum and position operators from the corresponding c-numbers) 

H C , 3D = p 2 /2m + hl{x) {|2)(l|e^ + h.c.} 

-^(2A + z 7 )|2)(2|, (1) 

where 7 is the Einstein coefficient of the excited level (level 2), i.e. its decay rate or 
inverse life time, Q(x) is the position dependent Rabi frequency (assumed to be real), 
A = uol — uj is the detuning, and p is the momentum operator in three dimensions (3D). 
The factor e %kh ^ takes into account the spatial dependence of the laser coupling. A one 
dimensional model is obtained by assuming that the atomic wave packet is centered at 
y = and satisfies k^Ay -C 1 before the first photon emission, so that the exponentials 
can be dropped and a one dimensional kinetic term suffices 0, 

.9 / h f n(x) \ 

where ft is the momentum operator conjugate to x, the ground state |1) is in vector- 
component notation m, and the excited state |2) is kJ. The probability, N t , of no 
photon detection from t , the instant when the packet is prepared far from the laser and 
with positive momenta, up to time t, is given by jlj 

N t = \\e- iH ^l h m ))\\\ (3) 

where \i/j(to)) is the (two-component) wave vector at t , and the probability density, 
Il(t), for the first photon detection by 

n(t) = ~ = 7^ 2 , (4) 
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where P 2 is the population of the excited state. 

To obtain the time development under H c of a wave packet incident from the left 
we solve first the stationary equation 

H C * = E<S>, where *(*) = (5) 

for scattering states with real energy E = h 2 k 2 /2m = E k , which are incident from the 
left (k > 0), 



/ e ikx +Rie -ikx\ 

\ R 2 e- i i x J 



x ~ — OO, 



V27T [ [t^) , X ~ OO. 

These states are not orthogonal, in spite of the reality of E, because the Hamiltonian 
H c is not Hermitian. The wavenumber q obeys 

E + = h 2 q 2 /2m, (7) 

with Img > for boundedness, while R\ t 2 and are reflection and transmission 
amplitudes for the ground and excited state channels. 

If ^j(k) denotes the wavenumber amplitude that the wave packet would have as a 
freely moving packet at t — 0, then 

V(x,t) = dki){k)^ k {x)e- lhkt/2m (8) 



n 

describes the "conditional" time development of the state for an undetected atom which 
in the remote past comes in from the left in the ground state. 

Some simple forms of Q(x) admit analytical solutions for the stationary waves, as 
demonstrated in [TJ |3] , but we shall limit the present discussion to an approximation 
that is valid for "arbitrary" shapes of Q within weak driving and low energy conditions 
0. 

Equation © reads explicitly 

tl 2 d 2 ^ 1) (x) s W9W x nI ,,W \ 

^ d J + ^(x)& 2 \x) = E$W(x), (9) 

i,2 02,^(2)/ \ * t 

+ -n( x )& 1 )( x )--(2A + t 1 )& 2 \x) = E& 2 \x). (10) 



2m dx 2 2 v ' y ' 2 
In the large 7 limit, 

' 2 Jl »-V,E, (11) 
we may neglect the kinetic and energy terms in the second equation to write 

g(2) = J 1 $(1) (12) 

2 A + z'7 v ; 

Physically, it is assumed that the excited state amplitude is small and proportional to 
the ground state amplitude, because the depletion by decay is rapid with respect to the 
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atomic motion and to the Rabi pumping period. Substituting (fT2|) into ©, there results 
a closed, one channel equation for $W with a complex potential 

= n 

v ; 22A + i 7 v ; 

~ 4A 2 + 7 2 ~ ' 4A 2 + 7 2 ' ( ^ 

so that the sign of the real part has the sign of the detuning whereas the imaginary part 
is negative, i.e. absorbing, for all x. 

The complex potential for the longitudinal direction x and a traveling laser wave 
perpendicular to the initial atomic motion is similar to the one that appears for the 
transversal direction y using the Raman-Nath or related approximations to describe 
the incidence of the atom on perpendicular standing laser fields EJ. In these 
approximations, however, x is treated as a parameter according to x = vt but y is 
kept as a variable, whereas here x is the variable and we consider the motion along a 
ray of constant y = 0. 

Notice that, within the stated conditions, the potential is independent of E. This 
implies that all (low energy) stationary scattering functions are subject to the same 
potential and therefore the time dependent Schrodinger equation also reduces to an 
effective one channel equation with the effective potential V. 

Within the one channel approximation, we have to find scattering solutions with 
the asymptotic form 



i [ Akx I d „—ikx 

$ (1) (*) = ^= t ' " " ~" (15) 

^2tt 1 T ie lkx , x~oo. 



The absorption (i.e. detection) probability A(k) for the incident &;-plane wave is given 
by 

A(k) = 1 — |i?i| 2 — | T\ | 2 . (16) 



3. Optimization 

The main aim of the present paper is to show that the atomic detection may be improved 
by varying the spatial dependence of the laser intensity and the detuning. To get 
optimal dependences for a fixed total length of the laser illuminated region and for a 
given momentum range one has to find first, for weak driving conditions, the form of 
the complex potential that maximizes absorption. 

A very similar objective (optimizing a complex absorbing potential) is also 
pursued in time-dependent molecular scattering calculations to eliminate the outgoing 
wave packets at the edge of the computational box and avoid unphysical effects 
E3 HU E21 nTn . We may take advantage of this coincidence by using, mutatis 
mutandis, similar optimization techniques. An obvious difference with the molecular 
scattering case is that here the real and imaginary parts of the potential have specific 
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physical meanings in terms of laser parameters. In particular the detuning A and Rabi 
frequency Q are given by 



This means that not all absorbing potentials V can be admitted, since the weak driving 
condition has to be satisfied, and moreover real and imaginary parts are not independent 
according to Eq. (|17|). The limitation of the weak driving condition may be avoided, 
but at the price of loosing the simplified one-channel description, i.e., by keeping the 
two-channel equations and a 2 x 2 complex potential matrix. In the same manner one 
could also get rid of the limitation to small kinetic energies and optimize detection of 
thermal atoms, for example. Here we shall discuss the simple one-channel case. 

The strategy is variational: we choose a functional form for V that depends on a few 
parameters and then find the values of these parameters that maximize the absorption 
for a a given velocity interval, or, in practice, for a discrete set of n velocities in such 
an interval [Zj, 



The "weights" W(kj) may be chosen according to the momentum distribution of the 
atoms, or, as in the examples discussed below, uniformly for a chosen absorption window, 
W(kj) = 1/n In principle, it is possible to construct explicitly potentials that absorb 
perfectly at a discrete set of wavenumbers [TT] . but these potentials tend to be too 
sensitive to small variations and thus ineffective for practical usage, and require in 
general arbitrary variations of real and imaginary parts. A further drawback in the 
present application is that they may have wild spatial variations or may correspond 
locally to strong driving conditions for which the one-channel, complex-potential model 
is not physically valid. 

Moreover, an ideal potential functional form for optimizing atom detection should 
enable us to control the spatial scale in which the laser intensity varies significantly, in 
accordance with technical capabilities. A simple form would be the sum of contiguous 
Gaussian functions with a certain width. Independently of the optimization algorithm 
used many evaluations of A are required in general. Each of them requires to solve 
numerically the Schrodinger equation to obtain the amplitudes T\ and R\, and this may 
be very time consuming. Here we shall choose an even simpler functional form, a set 
of contiguous "square barriers", as an approximation for a more realistic combination 
of Gaussians. At the present stage, where we are more interested in illustrating 
the concepts and general features involved than in any particular application, this 
simplification does not introduce any significant distortion The great numerical 
advantage of the square barriers is that the solution of the Schrodinger equation may 
be performed by multiplying a few 2x2 transfer matrices, which is an extremely 



A 



7 Re(V) 
2\m(V) ' 



(17) 






(19) 



Optimal atomic detection 



7 



fast process in comparison to the numerical techniques required for other functional 
forms. Moreover one can write explicitly the gradient of A with respect to the real and 
imaginary parts of the barriers, so that very efficient optimization algorithms may be 
used. Expressions of the absorption and its gradient in terms of transfer matrices may 
be found in |TTJ El EH] The main technical novelty here with respect to those works 
is the need to constrain the real and imaginary parts according to Eq. (fT7|) and to 
weak driving conditions. Our subroutine for constrained optimization is based in the 
successive quadratic programming algorithm [T4*j . 

4. Results and discussion 

Figure 121 shows the best absorption curves found for a range of momenta between 0.2 
cm/s and 9 cm/s considering one and two barriers for a total length of 10 [xm of the 
laser illuminated region, as in Figure 1. The potential, detuning and Rabi frequency 
corresponding to the optimal two-barrier case are shown in Figures E] and |3J Note that 
the optimal single barrier fails at both edges of the momentum range chosen whereas 
two barriers, with three free parameters to optimize, do a much better job. We have 
also calculated an optimal eight barrier potential, see figures El El and El that suggests 
that the form of an ideal continuous potential would be characterized by a smooth, 
non-linear increase of laser intensity with position, and negative detuning. 

The present variational approach is based on a number of simplifying assumptions 
but the concepts involved are applicable even outside the domain of validity of these 
assumptions. In particular, the restriction to weak laser interactions and low energies 
is not fundamental and may be removed by considering the full 2x2 potential matrix 
instead of the one-channel effective potential. Similarly, more realistic and smooth 
functional forms may be used. We have assumed here that a number of laser beams 
with different intensity may be combined in a composite double (or multiple) beam to 
maximize the detection probability, but there are other simple possibilities to explore: 
Meneghini et al., for example, have considered a linear modulation of the detuning, that 
could be realized in a nonhomogeneous magnetic field, and a standard Gaussian form 
for the Rabi frequency This amounts to a four-parameter functional form for the 
potential, whose absorption may be also maximized for specific applications. 

Apart from maximizing detection, there are other quantities that could be 
maximized or minimized, such as the detection delay. A minimal detection delay 
would be of interest for an accurate measurement of arrival times PU El ESI- One 
possible application of a "perfect detector" in a broad momentum range, would be 
the measurement, for the first time, of the backflow effect, namely, negative current 
densities for an atomic wave packet composed by positive momenta jU [TBI El CHI Ej- 
This and other extensions of the present work will be dealt with elsewhere. 
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Figure 2. Absorption versus velocity v for potentials optimized taking n = 100 
between 0.2 and 9 cm/s in i|19|) . L = 10 /im. Solid line: two barriers (each of 5 /mi); 
dashed line: one barrier; dots: eigth barriers. 
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Figure 3. Real (solid line) and imaginary part (dashed line) of the two-barrier 
potential of Figure 2. 
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Figure 4. Detuning (dashed line) and Rabi frequency (solid line) of the two-barrier 
potential of Figure 3. 
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Figure 5. Real (solid line) and imaginary part (dashed line) of an optimized eight- 
barrier potential. Other parameters as in figure |2 
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Figure 6. Detuning (dashed line) and Rabi frequency (solid line) of the eight-barrier 
potential of figure [3] 
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